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Abstract 

A novel approach to the problem of partial state estimation of nonlinear systems is pro¬ 
posed. The main idea is to translate the state estimation problem into one of estimation 
of constant, unknown parameters related to the systems initial conditions. The class of 
systems for which the method is applicable is identified via two assumptions related to the 
transformability of the system into a suitable cascaded form and our ability to estimate 
the unknown parameters. The first condition involves the solvability of a partial differen¬ 
tial equation while the second one requires some persistency of excitation-like conditions. 
The proposed observer is shown to be applicable to position estimation of a class of elec¬ 
tromechanical systems, for the reconstruction of the state of power converters and for speed 
observation of a class of mechanical systems. 
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1. Introduction 

The problem of designing observers for nonlinear systems has received a lot of attention 
due to its importance in practical applications, where some of the states may not be available 
for measurement. The interested reader is referred to for a recent review of the 

literature. 

In this paper a new framework for constructing globally convergent (reduced-order) 
observers for a well-defined class of nonlinear systems is presented. Instrumental to this de¬ 
velopment is to formulate the observer design problem as a problem of parameter estimation, 
which represents the initial conditions of the unknown part of the state. This new family 
of observers are called parameter estimation-based observers (PEBO). The class of systems 
for which PEBO is applicable is identified via two assumptions. The first one characterizes, 
via the solvability of a partial differential equation (PDE), systems for which there exists a 
partial change of coordinates that assigns a particular cascaded structure to the system that 
permits to obtain a classical regression form involving only measurable quantities and the 
unknown parameter. The second assumption pertains to our ability to consistently estimate 
this unknown parameter that, in general, may enter nonlinearly in the regression form. Eor 
linear regression forms, which may be obtained via over-parameterisation of the nonlinear 
regression, many well-established parameter estimation algorithms are available and the 


*This paper is submitted to Systems & Control Letters Journal. The abridged version of this paper will 
be presented at Conference on Descision and Control 2015, Osaka, Japan. 
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second assumption can be replaced by the well-known persistency of excitation (PE) con¬ 
dition |20l [27] . The latter condition is related with the “regular input” or “universal input” 
assumptions imposed in standard observer designs [a El do]. 

It should be underscored that, in contrast with the classical observer design method 
based on linearization up to output injection where the PDE to be solved imposes 
stringent conditions on the system, this is not the case for our PDE. The proposed PEBO 
also compares favourably with Kazantzis-Kravaris-Luenberger observers m in the follow¬ 
ing sense. Although both observers require an injectivity condition, in our observer this is 
imposed only on the partial change of coordinates mapping while in the Kazantzis-Kravaris- 
Luenberger observers the stronger requirement of injectivity of the full-state change of co¬ 
ordinates is needed. As is well-known [2] ensuring the latter injectivity property is the main 
stumbling block for the application of this kind of observers. 

The method is shown to be applicable for position estimation of a class of electromechan¬ 
ical systems. This class contains, as a particular case, the interesting example of permanent 
magnet synchronous motors (PMSM) that have been widely studied in the control and 
drives literature—see mm and references therein. It also allows us to design observers for 
a class of power converters under more realistic measurement assumptions than the existing 
results obtained with other observer design techniques. Finally, it generates simple speed 
observers for mechanical systems that are partially linearisable via change of coordinates 
(PLvCC)—a practically important class that has been thoroughly studied in [2^ . 

The remaining of the paper is organized as follows. Section presents the problem 
formulation and main result. Section El is devoted to a discussion of the results. The case 
of linear time-invariant (LTI) systems is treated in Section Section illustrates the 
application of the technique to three physical examples. The paper is wrapped-up with 
concluding remarks and future research directions in Section]^ 


2. Problem Formulation and Main Result 


In this section the (partial state) observer problem addressed in the paper, and the 
approach that we propose to solve it, are presented. The class of systems for which the 
PEBO design technique is applicable is identified via two assumptions. The first one. 


given in Subsection 2.2, characterizes systems for which there exists a partial change of 
coordinates that assigns a particular cascaded structure to the system that permits to 
reformulate the state observation problem as a problem of parameter estimation. The 
second assumption, given in Subsection |2.3[ pertains to our ability to consistently estimate 
this unknown parameter. 


2.1. Partial state observer design problem 
Consider the dynamical system 


X = fx{x,y,u) 

y = fy{x,y,u), (I) 

where fx ■ x x —>■ and fy : x M”*' x —)■ M"'?' are smooth mappings 0 

Assume that the input signal vector u : M+ —>■ is such that all trajectories of the 

system are bounded. Find, if possible, mappings F : x x —>■ and G : 


^Throughout the paper it is assumed that all mappings are sufficiently smooth. 
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X X M™' —)• , for some positive integer , such that the (partial state) observer 

? = F{^,y,u) 

x = G{C,y,u), ( 2 ) 

ensures that ^ is bounded and 

lim |x(t) — x{t)\ = 0, (3) 

t^OO 

for all initial conditions (x(0),y(0),.^(0)) G +"•!/+»*? and a well defined class of input 
signals u gU. 

It is important to underscore that, in contrast with the usual observer problem formu¬ 
lation, a provision regarding the input signal is added. This additional qualifier is needed 
because the observation problem will be recast in terms of parameter estimation whose solu¬ 
tion requires “sufficiently exciting” signals. See Section 1.2 in [5] for a thorough discussion 
of the role of the input in the observation problem. Also, note that we have writen the 
system dynamics including the output y as part of the state, this is done, of course, without 
loss of generality. 


2.2. System re-parametrization 
Assumption 1. There exists three mappings 

(/> : M"" X ^ M"" 

h:R^y xR^ ^ 

with riz > nx, verifying the following conditions. 

(i) (Left invertibility of •) with respect to its first argument) 

(l>^{(l){x,y),y) = X, Vx G , Vy G . 


(ii) (Transformability into cascade form) 


dcj) 

dx 


fx{x,y,u) + 


dcj) 

dy 


fy{x,y,u) = h{y,u). 


( 4 ) 

□□□ 


An immediate corollary of (ii) in Assumption is that the partial change of coordinates 


z = (l){x,y), 


( 5 ) 


ensures 

z = h{y,u). (6) 

Moreover, the left invertibility condition (i) ensures that the partial state x can be recovered 
from z and y, that is, 

x = 4>^{z,y). (7) 

The cascade structure of the system is given in Fig. 
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Figure 1: Block diagram representation of of the transformed system. 


Proposition 1. Consider the system Q verifying Assumption Define the dynamic 
extension 

X = h{y,u), (8) 

with x(0) G We can compute a mapping x x x —>■ M”'*' such that 

y = (9) 

X = (l)^{x + 0,y), ( 10 ) 

where 0 G is a vector of constant, unknown parameters. 

Proof 1. From ([^ and Q we get i = x- Hence, integrating this equation yields 

zit) = xit) + ^, ( 11 ) 

where 


9:= ^(O)-x(O). 


( 12 ) 


m 


Replacing (0 
in 0 to get 


Q yields (10). 


Finally, the regression model ([^ is obtained replacing 


fyi4>^ix + y),y, u) =■ ^(x, 2 /, u, 9), 

completing the proof. 


(13) 

□□□ 


2.3. Consistent parameter estimation 

An immediate consequence of Proposition is that the problem of observation of the 
unmeasurable state x is translated into a standard parameter estimation problem for the 
regression model 0 with the observed state generated by 

x = (t>^{x + 0,y), (14) 

where 9 : M_|_ —)■ is an on-line estimate of the constant vector 9. Therefore, to complete 

the PEBO design it is necessary to ensure the existence of a consistent estimator for the 
unknown parameter 9. Towards this end, the assumption below is introduced. 

Assumption 2. There exists two mappings 

H X X M”!' x M”" 

N X X R^y X M'" ^ 


with > 0 such that the parameter estimator 

C = H{x,C,y,u) 

9 = N{x,C,y,u), (15) 

coupled with the dynamic extension ([^ and the regression model @ ensures that C is 
bounded and 


lim 

t^OO 


9{t) - 9 


= 0 , 


(16) 


for all initial conditions (y(0), x(0), C(0)) £ and a well defined class of input 

signals u^lA. □□□ 
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2.4- Main result 


The main result of the paper is contained in the corollary below whose proof follows 
immediately from ([8|, (10), (14), (|15|) and the parameter convergence assumption (16). 


Corollary 1. Consider the system 0 verifying Assumptions and with u £ U. A 
(partial) state PEBO of the form ([^ that guarantees ([^ is given by 


Fi^,y,u) := 
G{C,y,u) := 


col(x,C) 

h{y,u) 

_ H{xX:y,u) 

4>^{x + N{x,C,y,u),y). 


(17) 

□□□ 


2.5. On the solvability of the PDE (Q 

Similarly to all constructive observer design methods the proposed technique involves 
the solution of a parameterised PDE, namely (Q— where we recall h{y, u) is a free function. 
See O m [5l [To] for a recent review of the literature where the PDEs of the various existing 
observer design methods may be found. A key step in our observer design is, of course, the 
explicit solution of this PDE, some comments in this respect are made in this subsection. 

To streamline the presentation of the subsequent discussion it is convenient to define 
the state vector X := col(x, y) and the (u-parameterized) vector fields 


UX) := eo\{UX,u)Jy{X,u)) 
hu{X) := h{y,u), 

and rewrite the PDE (Q element-by-element as 

dUX) 


dX 


-fu{X) = {hu)i{X), i = l,...,n^ 


Eollowing the construction used in [8| we then define Uz functions 

* — 1) • • ■ ) 


pnj;+n. 


V X 


(18) 


as 


and write the non-homogeneous PDE (18) as the following homogeneous PDE 

d^i{X,s) 


dX 


'fu{x) + ^^^^i^{huUx) = o. 


Erobenius Theorem m states that, if 

rank 


MX) 

{hu)iiX) 


ds 


= 1 , i = l,...,nz, 


(19) 


( 20 ) 


(uniformly in u) in a neighbourhood of Xq G then (19) has a loeal solution around 

a point {Xq, sq) G x M, because all one-dimensional distributions are involutive. 

Unfortunately, the regularity condition (20) may be a restrictive assumption in the present 
scenario, since rules out solutions around equilibrium points of the system ©• This renders 
Erobenius Theorem unaplicable for our problem in—often encountered—stabilisation tasks. 
See Remark R5 in the next section. 

Even in the case where the rank condition is satisfied it is clear that existence of so¬ 
lutions does not imply that an analytic expression for it can be easily obtained—and the 
requirement of solving the PDE remains the main stumbling block for our approach. 
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3. Discussion 


R1 Besides the explicit solution of the PDE Q an additional difficulty is the selection of 
a mapping h{y, u) that will ensure that the partial change of coordinates 4>{x, y) admits 
a left inverse (with respect to x). It should be noted that this injectivity-like property is 
imposed on the partial change of coordinates 4>{x,y), which should be contrasted with the 
requirement of injectivity of the full-state change of coordinates imposed in the Kazantzis- 
Kravaris-Luenberger observers m- As discussed in |2] the latter injectivity property is 
the main stumbling block for the design of these observers—see [26] for a very illustrative 
example. 


R2 For general nonlinear systems the regression system Q depends nonlinearly on the 
unknown parameters. Although some results are available for the estimation of nonlinearly 
parameterised, nonlinear systems o m na da |28] the problem of generating consistent 
estimates remains wide open. On the other hand, for the case of linear parameterisation 
the estimation problem has a standard solution. Indeed, many techniques lacniEii are 
available to generate consistent estimates for linear regressions of the form 

y = ^o{x,y,u) + ^i{x,y,u)e, ( 21 ) 

with known mappings <i>o : x M”?' x —)■ R^y and <hi : x x —)■ 

See the examples in Section 


R3 As is well-known |20|, the parameter convergence requirement in parameter estimators 
involves some form of excitation on the signals—this requirement is encrypted in the condi¬ 
tion u €U of Assumption This condition is of the same nature as the “universal input” 
or “regular input” conditions for classical observer designs HE] do). For linear regressions 
of the form ( |21[ ) it has a very precise characterisation in terms of PF of the regressor matrix 
<hi(x, y, u), which is defined as the existence of constants d > 0 and T > 0 such that for all 


t > 0 


rt+T 


^Jixis),y{s),u{s))<f>i{xis),yis),u{s))ds > 5In 


( 22 ) 


Under the PF condition above it is straightforward to design globally convergent parameter 
estimators for the regression model (21). In this case the set lA is defined as follows: 


U ■.= {u : M+ 


(22) holds along trajectories of (j^, ([^}. 


From (22) it is clear that if there are more measured states than unknown ones, that is, 
if Uy > Uz, then the PF condition translates into a simple rank condition on the matrix 
<hi(x, y,tt)—this is the case of mechanical systems treated in Subsection 5.1 


also recalled that in the identification literature there are well-known relationships between 
the adaptation gains, the PF constants 5 and T and the convergence rate of the estimation 
errors; see iniEii. 


R4 It should be underscored that, in many cases, it is possible to transform a nonlinearly 
parameterised regression into a linear one via over-parameterisation —see the discussion in 
this respect in m- Since over-parameterisation increases the dimension of the parameter 
space the excitation requirements on the signals are, of course, more stringent. 


R5 A potential practical drawback of the proposed technique is the utilisation of pure 
integrators in ([^. Indeed, in some applications the measurable signals may exhibit a (sign 
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definite) bias in steady-state that will lead to unbounded signals when fed into open-loop 
integrators. On the other hand, this problem is conspicuous by its absence in regulation 
tasks. Indeed, in this case there exists a desired, constant operating point (x*, y*) G 
that must satisfy the equilibrium equations 


0 = fx{x*,y*,u^) 

0 — 

for some constant tt* G M™. From the equations above it is clear that a necessary condition 
for solvability of the PDF Q is that h{y^,u^,) = 0. Hence, in normal operating conditions, 
the open integration operation will not generate a bias. A similar scenario appears in the 
ubiquitous PI controllers widely used in industry to drive some error signal to zero. To 
avoid drift— e.g., in the presence of noise—several ad hoc remedies, including the addition 
of small leakages and resettings, are well established. 


R6 As shown in (12) the nnknown parameter 0 —and, consequently, the estimated state—is 
determined by the system and observer initial conditions. It may be then argued that this 
makes our analysis “trajectory dependent”, hence intrinsically fragile. This criticism would 
certainly be pertinent if off-line, instead of on-line, parameter estimators were advocated 
or if transient performance claims were made. Since this is not the case in the present work 
the argnment seems specious. 


4. Case of Linear Time—Invariant Systems 

The proposed observer design procedure is of little—if at all—use for LTI systems. 
However, it is interesting to show that even for this simplest case the relationships between 
the classical notions of observability [Hj and identifiability [30] and Assumptions and 
are far from obvious. 


Proposition 2. Assume the system ([^ is LTI, that is. 


X 


All 

Ai 2 


X 


'Bi 


y. 

— 

A21 

A22 


y. 


.^2. 

u 


(23) 


where Aij and Bi, i,j = 1, 2, are constant matrices of suitable dimensions. 


(Cl) Observability of the system (23) (with respect to the output y) does not imply As¬ 
sumption 

(C2) Assumption!^ does not imply observability of the system (23). 

(C3) If Assumption holds then observability of the system (23) is necessary for identifi¬ 
ability [30] of the parameter 9 (defined in Proposition [T). Hence, it is necessary for 
Assnmption to hold. 


Proof 2. Without loss of generality we take 


4>{x,y) =Tix + T 2 y, 


where Ti G and T 2 G Since Uz > Ux, 

ensured imposing 

rank Ti = Ux- 


condition (i) of Assnmption is 

(24) 
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Under this conditions we have 


4>^{z,y)=Tl{z-T2y) 


(25) 


where 

tI := {T'JTi)-^T'^. 

The PDE (|^ reduces to the algebraic equation 

Ti^ii + T 2 A 21 = 0 . (26) 

Now, applying Popov-Belevitch-Hautus criterion m we conclude that the system ( |23| ) 
is observable if and only if 


rank 


All 

A2I 


— 


Vs E a{Aii}, 


(27) 


where fT{-} denotes the set of eigenvalues. 

We prove claim (Cl) by constructing a system that is observable but does not satisfy 
Assumption!^ For, take Ux = 2 and Uy = 1 and set An = diag{ai,a 2 } with 01 / 02/6 
and A 21 = [ a ^ ] with o, 6 / 0. The matrix appearing in (27) takes the form 


sh - An 


s — Ol 

0 

A 21 

— 

0 

0 — 02 
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a 

b 


whose rank is two for s = Oj, i = 1, 2. Hence, the system is observable. On the other hand, 


for all > 2 , (26) is equivalent to 


Ti — —T2A2 iA^^ . 


Consequently 


rank Ti < minjrank T 2 , rank A 21 } = 1 < 2 = Ua: 


violating the rank condition (24) 


Claim (C2) is proven by contradiction constructing a system that is not observable but 
satisfies Assumption For, take Ux = Uy = Uz = 1 and set A 21 = 0 that—from (27)— 


implies the system is not observable. However, if An = 0 the algebraic equation (26) admits 
a solution Ti = 1 , T 2 = 0 ensuring Assumption [Tj 

Finally, claim (C3) is established proving, after some lengthy but straightforward cal¬ 
culations, that the regression form ([^ is given by 

ij = A 21 TIX + (A 22 - A 21 TIT2)y + B 2 U + A 21 TI0, 


where we notice that A 2 ir/ E while A 21 E 

identifiable if and only if Uy > Uz and rank A 21 = Ux — 
it is clear that observability follows. 


From m it follows that 6 i 


IS 


from the latter condition and (27) 


5. Application to Three Physical Examples 

In this section we prove that PFBO is applicable to the speed observation of PLvCC 
mechanical systems studied in |29j , the position observation for a class of electromechanical 
systems, and the reconstruction of the full state from partial measurements of a popular 
switched power converter. 














5.1. Mechanical systems which are partially linearisable via change of coordinates 

In this subsection we are interested in the problem of speed observation of mechanical 
systems described in Hamiltonian form by 


y 

X 



r dH -\ 




dy 


0 


dH 

+ 

. G{y) 


- dx - 




(28) 


where s \= ^ is the number of degrees of freedom of the system, y, x G M* 
position and momenta, respectively, u G is the control input, m < 
is the full rank input matrix. The Hamiltonian function H x 
function 

1 


H{y,x) = ^x'^M ^{y)x + V{y), 


is the positive definite inertia matrix and V : 


are the generalised 
s, G :W ^ 

—)> M is the energy 


is the potential 


where M : M® —)> S 

energy function. It is assumed that position y is measurable and we want to estimate 
velocity y via the estimation of momenta x and the relation y = M~^{y)x. 

It will be shown that if the dynamics can be rendered linear in momenta (velocities) via 
a change of coordinates then Assumptions 1 and 2 of the paper are satisfied. Moreover, the 
resulting reparameterisation is linear, that is, of the form (21), and the PE condition (22) 
is trivially satisfied, hence the set Id is the whole input space. These systems, referred as 
PLvCC, have been studied in |29] and they have been characterised via the solvability of a 
PDE—see also [7] for an intrinsic characterisation of the class. 

To present our result we need to recall the following assumption from [29]. 


Assumption 3. Given the inertia matrix M{y). There exists a full rank matrix T : M® —)> 
M®^® such that, for f = 1,..., s. 


where the matrices : 


H(,)(y)+Hj)(y) = 0, 
^sx® defined as 

1 






(29) 


(30) 


with Ti{y) the f-th column of T{y), Tij{y) its ij-ih. element and \Ti,Tj\ the standard Lie 
bracketIn this case the mechanical system (28) is PLvCC. 


Proposition 3. Consider the mechanical system (28) whose inertia matrix verifies As¬ 
sumption 

(i) Assumption is satisfied with the mappings 


fix.y) = {y)x 

(t>^{z,y) = T^"'{y)z 

h{y,u) = -T^iy) 


dV 

^{y)-G{y)u 


standard Lie Bracket of two vector fields Ti{y), Tj{y) is defined as [71,7^] ;= ^^71 — ^Tj. 
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( 31 ) 


(ii) The mapping <h(x, y, of Propositionis linear in 6 and yields 

y=[T^{y)M{y)]-\x + 0)- 

(iii) The PE condition (22) is satisfied for all input signals u. Hence, Assumption holds 
with U being the whole input space. 

Proof 3. Invoking Proposition 1 of [221 have that the partial change of coordinates 

= T'^{y)x, 


transforms the system into the form 

y = [T^ iy)M{y)]~^z 
z = -T^{y) 


dV 

^{y)-G{y)u 


which establishes claims (i) and (ii). To prove claim (iii) we refer to (21) and (31) and 
identify 

<^i{x,y) ■■= [T^{y)M{y)]-^ 


which is a square, full rank matrix, hence (22) is trivially satisfied. 


From Proposition the design of a globally exponentially stable momenta observer 
follows trivially, and is omitted for brevity. It should be underscored that the resulting 
observer is much simpler than the I&I observer proposed in |29j . This is particularly true 
for systems that do not satisfy the integrability Assumption 2 in [29] , for which it is necessary 
to use the, rather involved, and high-gain-\ike technique of I&I with dynamic scaling. 


5.2. Position observation in electromechanical systems 

We consider electromechanical system^ consisting of nx inductances and a single mass, 
whose flux and mechanical position are denoted by A G and q G M, respectively. The 
magnetic energy stored in the inductances is given by 

£n{i,q) ■■= ^i^L{q)i + {q)i 

where L : M —)• is the positive definite, (position-dependent) inductance matrix, 

/i : M —)> represents the flux linkages due to the possible presence of permanent magnets 

and i G are the currents flowing through the inductances. The kinetic energy of the 

mass is 

i^{q) ■= 

with j > 0 the (constant) mass inertia. We assume that the system is subject to constant 
external forces, e.g., gravitational forces, but it does not have any other potential energy 
storing elements. Hence, the systems potential energy is given by 

V{q) := qr, 

where r G M represents the constant external force. 


^The interested reader is referred to | 22l I24| for additional details on this model. 
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The dynamical model of the system is obtained applying Euler-Lagrange equations with 
the energy functions above yielding 


L{q)—+L'{q)qi + ^'{q)q + Ri = Bu 
jq-^i^L'{q)i-i^fi'{q) +fq = -r, 


(32) 

(33) 


where (•)' denotes differentiation, R = diagjri,..., Vn^} > 0 is the matrix of resistors (in 

are external 


series with the inductors), B G constant input matrix, u G 

voltage sources and / > 0 is a Coulomb friction coefficient. 

As shown in 


the model (32), (33) describes the behaviour of a large class of elec¬ 
tromechanical systems, including the classical levitated ball and the most common electrical 
motors. 

We are interested in the design of an observer for the mechanical position q measuring 
col(z, q), as well as the more practically interesting case when we measure only i. To proceed 
with the observer design we recall Gauss’s and Ampere’s laws that establish the following 
expression for flux linkage vector 


A = L{q)i + ^i{q). 


(34) 


Moreover, Gauss’s law tells us that 


A = —Ri + Bu. 


(35) 


The latter two equations correspond, of course, to the electrical equation (32). Since i 
is measurable ( [3^ proves that L{q)i + ^{q) qualifies as an admissible partial change of 
coordinate^ (^(x, y) verifying condition (ii) of Assumption with Uz = 2 and 

h{y, u) := —Ri + Bu. 


Gondition (i) of Assumption is satisfied if—given (34)—we can recover q from measure¬ 
ment of A and i. Interestingly, we show below that this is the case for PMSMs. 

Unfortunately, if the only measurable quantity is i, its time derivative ^ in (32) is not 
in the form of (Q, that is y = fy{x, y, u), with x being only q and y only i. Indeed, besides 
q and i, the electrical equation (32) contains the velocity q. Hence, in order to obtain the 


regressor from ([^ it is necessary to assume also measurement of q. 

Let us proceed now with the observer design for the (surface mount) PMSM |16l I24| . 
For the sake of clarity of exposition assume first that i and q are measured—the requirement 
of measuring q is relaxed later. For the PMSM we have nx = 2, m = 2, B = I 2 and 


m 


Lh 

Am 


cos(npg) 

sin(npg) 


(36) 


where the positive constants L, Am and Up are the stator inductance, permanent magnet 
flux constant and number of pole pairs, respectively. Hence, defining 


(t){q,i) 


Li -|- Am 


cos(npg) 

sin(npg) 


A, 


^Notice that, with respect to the notation in PropositionA plays the role of z and q is x. 
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it is clear that the mapping 


:= —arctan 

Up 

satishes condition (i) of Assumption 

The dynamic extension ® is given by 


A 2 — Li2 

Ai — Lii 


X = —Ri + u. 


Now, from (35) (with B = I 2 ) and (37) we have, upon integration, that 

Kt) = x{t) + 0, 


( 37 ) 


(38) 


where 0 := A(0) — x(0) is the unknown parameter. After some lengthy, but straightforward 
calculations mimicking the proof of Proposition we obtain a linear regression form for 
the currents as 


dt 


= + $i(g)6» 


(39) 


where 


‘^oix,i,Q,u) := 


L 

UpQ 


-UpQ 

L 


+ Upq 


X2 

Xi 


+ u 


$i(g) := UpQ 


X2 

XI 


Although the regression form (39) can be extended with an equation for q this turns out to 
be unnecessary to solve the parameter estimation task that consists only of two unknown 
parameters. 

A classical parameter estimator can be designed for the linear regression form (39). 
However, we make the important observation that the requirement of measuring q, which 
is not realistic in a practical scenario, can be obviated. Towards this end we proceed as 
follow. First, from (34) and (36) we have that 

|A-Lf = A^. 


Second, replacing (38) above and expanding the square yields the static linear regression 
form 

Yix,'!-) = (40) 

where 

y{x,i) ■= \X-Li\‘^ 

- -2{x-Li) ^ 


S{x,i) ■= 


1 


are, of course, measurable and the new (extended) unknown parameter is 


rj : = 


9 

- l«P 


A full theoretical analysis and extensive simulations and experimental results of parameter 
estimators for the regressions (39) and (a filtered version of) (40) may be found in [6]. As 
shown in that paper the set lA is defined as follows: 

rt-\-T 

U ■.= {u ■. M+ —I / q^{s)ds > (5 > 0, along trajectories of (32), (33) and (36)}. 
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5.3. Cuk converter 

In this subsection we apply the proposed observer design technique to power converters. 
As will become clear below the technique applies to a broad class of converters, including 
the popular boost converter. For the sake of ease of exposition, instead of developing 
a—notationally cumbersome—general theory for a broader class of power converters, we 
preferred to concentrate on the specific example of the Cuk power converter, depicted in 
Fig.m 


Li 6^2 



The average model of this device is given by 

Li<j^ = -{l-u)v 2 + E 
C2V2 = - u)ii + uiz 

^3^ = -UU2-U4 ^ ^ 

C/^V4, = is — Gv4, 

where Li, ( 72 , T 3 , C 4 , E and G are positive constants and u G (0,1) is a duty cycle. We refer 
the reader to [ 1 ] for further details on the model. 

To illustrate the generality of the approach we consider two different measurement sce¬ 
narios. In the first one we assume that (^ 2 ,^ 4 ) are measurable, while in the second one 
( 02 , 13 ) are measurable. Although from the practical viewpoint it is “easier” to measure 
voltages, we also consider the second one since, as shown in [ 1 ], is the one that can be 
solved with immersion and invariance (I&I) observers, with which we compare our observer 
in simulations below. 


Case I Denoting y := 


col(u 2 ,U 4 ), X := col(ii,i 3 ) we get from (41) 


Xl = 

X 2 = 


- —(1 - u)yi + 


L 3 


uyi 


L 3 


2 / 2 - 


E 

Ti 


Since the right hand side of these equations is independent of x we can directly select 

4>{x,y) = X. 


The dynamic extension is given by 


X = 


E 1 


\uyi - ^^y2 


La' 


=: h{y,u), 


and the regression form is 


y = ^o{x,y,u) + ^i{u)e 

e := x( 0 )-x( 0 ) 


( 42 ) 
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where 


^oix, y,u) ■= 


'"^(1 + d ^ uX 2 

- gy2 


, ^>i(u) : = 




0 


C4 


(43) 


The model (42) contains the time derivative of the output y. To get a classical (static) 
regression model we use the standard filtering technique |23j and define the filtered signals 


a 


p + a 


(•), 


where P '■= and a > 0 is a design parameter. Applying the filter to (42) we obtain the 
standard linear, static regression model 


'd = + e 


(44) 


where 


:= 


ap 


a 


-y 


-4>n 


p + a p + a 

is clearly measurable (without differentiation) and e is an exponentially decaying signal that 
depends on the filter initial conditions and the filter time constant 


The regression model (44) is used for the parameter estimator, which is the classical 
gradient estimator 

e = r¥7(79-^i0), 

where the adaptation gain T = T^ > 0 is a design parameter. The state observer is defined 

as X = 9 + X- 


It is important to underscore that the regressor matrix 4>i(u) given in (43) has an 
extremely simple form. Indeed, due to its upper triangular form, the estimation of the 
second parameter is decoupled from the first one and, moreover, the corresponding term in 
the regression is simply the constant Also, since the matrix depends only on the input 
signal u the set U is defined as 


ZY := {u : M+ ^ (0,1) 


rt+T 


1 — u{s) 
{l-u{s))u{s) 


{l-u{s))u{s) 


n^) + 


ds > 6 I 2 > 0 }. 


Some simple calculations show that the matrix inside the integral is positive definite for any 
u G (0,1). Hence, U = {u ■. M+ —)■ (0,1)} and consistent estimation is always guaranteed. 
Simulations were carried out to evaluate the performance of the proposed observer. The 


simulations were done for the model (41) in closed-loop with the certainty equivalent version 


of the full-state feedback I&I controller given in Proposition 8.2 of [1]. That is, the control 
law was defined by 


u = 




+ A: 


G\Vd\v 2 + E{x2 - Xi] 


\Vd\+E 1 + (G\Vd\v2 + E{x 2 - xiW 

where < 0 is the reference imposed to the output voltage U 4 and A is chosen as 


(45) 


A = Ao min 




E 


\Vd\+EfiVd\+E 


with 0 < Ao < 2. The full-state version of this controller, i.e., replacing xi and £2 by ii 
and Za, respectively, ensures global asymptotic stability of the desired equilibrium as well 
as verification of the saturation constraints in the input signal. 
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The numerical simulations were performed with the following values of the converter 
parameters Li = 10 mH, C 2 = 22.0 ;uF, L 3 = 10 mH and C/^ = 22.9 /rF, G = 0.0447 S 
and E = 12 \. The initial conditions for all simulations are set to x(0) = (0.5,—1), 
y(0) = (10,-12). The initial set point for the output voltage is Vd = 25 V, and then 
this is changed at t = 0.2 s to = 30 V, at t = 0.4 s to = 15 V, at t = 0.6 s to = 5 V, 
at t = 0.8 s to Vrf = 20 V. The simulation results are presented in Fig. 


Case II Denoting now y := col(u 2 ,i 3 ), x := col(ii,U 4 ) we get from 


1 , , E 

xi = -^(1 - u)yi + ^ 

Li Li 


1 


G 


X2 = 7r2/2 - ^X2. 

04 04 

The right hand side of the second equation depends on X 2 , therefore the choice 4>ix,y) = x 
is not suitable here. We propose instead 


0 


4>{x,y) = X - 

that, introducing the partial change of coordinates z = 4 >{x,y), yields the required form 

= : h{y,u). 


z = 


^2/2 + guyi 


The dynamic extension is then given hy x = h{y,u), and the regression model is of the form 

y = ‘^oix,y,u) + <^i{u)e 


e ■= a;(0)-x(0)- 

where 

^o{x,y,u) := 

The state observer takes the form 


0 

^ 2 / 2 ( 0 ) 


L C 4 


^(1 - u)xi + ^uy 2 
- 7^x2 - ^y 2 




4>i(u) : = 


0 

_ j_ 
Ls 


(46) 


X = 9 + X + 


0 

GLs 

Ci 


y2 


The regressor matrix 4>i(u) given in (46) has an even simpler form than the one of Case 
I above. Indeed, the matrix is now diagonal with the second term in the regression simply 
the constant Clearly, for this case we also have U = {u : M+ —)• (0,1)} and consistent 
estimation is always guaranteed. 

In Proposition 8.3 of [1] the following I&I observer is proposed 


xmi 

L 

C 2 


Cl 

C2 


C27iyi 

^372^2 


1 

1 


C 4 


[-(1 - u)yi + E]- 7 i[(l - u)(Ci + ^ 27171 ) + uy 2 ] 
[y2 - g{C2 - ^3722/2)] - 72 + C2 - ^3722/2], 


(47) 
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where 71,72 > 0 are design parameters. It should be noted that in the latter reference the 
parameters E and G are treated as unknown and are also estimated. If they are assumed 
known the I&I observer takes the form given above. 

The performance of our observer was compared with the I&I observer (47) via numerical 
simulations. They were done under the same scenario as the ones done for Case I, but 
now with the certainty equivalent observer that results replacing ii and V 4 by xi and X 2 , 
respectively. The simulation results are presented in Figs, [^[^with different observer gains. 


6. Conclusions 


A radically new approach to design state observers for nonlinear systems has been 
proposed. The key idea is to translate the state observation problem into one of parameter 
estimation. It turns out that this is possible if we can find a partial change of coordinates 
z = 4>{x, y) such that i depends only on y and u. The observer then comprises a copy of 
i (called x) that, upon integration, differs from z only on the initial conditions—and these 
are the (constant) parameters that we propose to estimate. If the change of coordinates 
satisfies an injectivity property then x can be estimated from the knowledge of y, u, x and 
an estimate of 9. Clearly, if the latter converges to 6 , then the estimate of x will converge 
to its true value. 

It has been shown in the paper that the change of coordinates is obtained from the 
solution of a parameterised PDF, which does not impose (a priori) the strict constraints 
of the classical observer design HZ! Moreover, it is argued that the required injectivity 
property is weaker than the one required in the Kazantzis-Kravaris-Luenberger and the 
I&I observers. 

The design of the observer is completed adding a parameter estimator to a regression 
model of the form ij = ^(x,?/,^, 0) that, in general, depends nonlinearly on the parame¬ 
ters. Although some estimation techniques for nonlinearly parameterised nonlinear systems 
are available, it is also suggested that—via over-parameterisation—it may be possible to 
transform the regression into a linearly parameterised one. The latter case has been widely 
studied in the literature and many algorithms that guarantee parameter convergence under 
some excitation conditions are available. 

The proposed technique has been shown to be applicable to position estimation of a class 
of electromechanical systems, to power converters and to speed estimation of the PLvCC 
mechanical systems studied in [29] . 

Current research is under way in the following directions. 


Identify other classes of physical systems to which the proposed method is applicable. 

Compare the performance of PEBO for the PLvCC mechanical systems discussed in 
Subsection O with other existing speed observers. 

Further clarify under which conditions the required partial change of coordinates exists 
and when it will lead to an easily tractable, linearly parameterised, regression model. 


• Further explore the connection between the classical concepts of observability and 
identihability and Assumptionj^and the excitation conditions required by the method. 
In this respect, the analysis of the simplest LTI case of Section [^ shows that this task 
is far from obvious. 
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(a) (b) 




(c) (d) 


Figure 3: Transients of the observation errors (a) xi := xi — ii , (b) X 2 '■= X 2 — * 3 , (c) the 
voltage reference Vd and voltage output U 4 and (d) the control input u for the tuning gains 
a = 0.5, r = 0 . 001 / 2 . 
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(a) (b) 




t, s 

(c) 


t, S 

(d) 


Figure 4: Transients of the observation errors (a) xi := xi — ii , (b) X 2 '■= X 2 — V 4 , (c) the 
voltage reference Vd and voltage output V 4 and (d) the control input u for the proposed 
observer with tuning gains a = 1, F = diagjO.Ol, 0.1}. 
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(a) 



(c) 



(b) 



(d) 


Figure 5: Transients of the observation errors (a) xi := xi — ii , (b) X 2 '■= X 2 — V 4 , (c) the 
voltage reference Vd and voltage output V 4 and (d) the control input u for the proposed 
observer with tuning gains a = 1, F = diagjl, 10}. 
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(c) (d) 


Figure 6 : Transients of the observation errors (a) xi , (b) X 2 , (c) the voltage reference Vd 
and voltage output U 4 and (d) the control input u for the I&I observer (47) with tuning 


gains 71 = 25, 72 = 1. 
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